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The  time  will  come  when 
diligent  research  over  long  periods 
will  bring  to  light  things  which  now 
lie  hidden.  A  single  lifetime  ... 
would  not  be  enough  for  the  investi¬ 
gation  of  so  vast  a  subject  ...  Many 
discoveries  are  reserved  for  ages  yet 
to  come  . . . 

Seneca,  Natural  Questions 
6ocm  7,  first  century 

As  is  so  often  the  case  in  scientific  research,  the 
experimenter,  in  his  desire  < :o  uncover  the  truth,  unearths 
more  questions  than  he  answer's.  I  can  now  empathize  with 
that  scientist  after  havl'  o  worked  for  so  long  on  this 
thesis.  As  data  was  slcvU  acquired,  so  too  were  theories  to 
explain  the  results,  Htwtver,  further  data  reduction  would 
often  prove  the  hypothesis  in  error,  and  a  re-thinking  was  in 
order.  Yet  when  data  and  theory  finally  agree,  you  enjoy  a 
wonderful  sense  of  accomplishment.  As  Democritus  of  Abdera 
once  said,  "I  would  rather  understand  one  cause  than  be  King 
of  Persia." 

None  of  this,  however,  would  have  been  possible  without 
the  help  of  several  individuals.  First  and  foremost  is  Dr. 
James  Hitchcock,  my  thesis  advisor.  I  first  became 
interested  in  heat  transfer  during  an  introductory  course 
which  he  had  taught.  Throughout  my  entire  heat  transfer 
curriculum,  I  have  always  been  eager  to  learn  from  and  study 
under  Dr.  Hitchcock.  Secondly,  I  would  like  to  thank  2Lt . 
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James  Shoemaker,  a  member  of  the  AFWAL  Propulsion  Lab,  for 
the  use  of  his  digitization  equipment.  Thanks  also  go  out  to 
Mr.  John  Brohas  and  the  AFIT  lab  technicians  for  their  help 
in  constructing  my  apparatus. 
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The  two  classical  problems  ot  charging  a  tank  with  an 
Ideal  gas  concern  themselves  only  wltn  the  initial  and  final 
states  of  the  gas.  The  first  problem  involves  an  insulated 
tank.  From  initial  conditions,  the  final  temperature  of  the 
gas  in  the  tank  can  be  determined.  The  second  problem 
involves  a  tank  with  conducting  walls.  In  this  case,  it  Is 
the  total  heat  transferred  from  the  gas  to  the  tank  walls 
which  can  be  found.  This  study  went  beyond  these  equilibrium 
analyses  to  acquire  a  fundamental  understanding  of  the 
transient  phenomenon  involved  In  charging  a  tank. 

By  way  of  experimentation,  transient  pressure  and  heat 
flux  measurements  were  made  on  an  evacuated  cylindrical  tank 
as  1 t  was  charged  by  sudden  exposure  to  the  atmosphere . 

Then,  an  average  convection  coefficient  and  driving 
temperature  potential  for  the  tank  were  found.  A  method  was 
then  developed  for  estimating  heating  rates  for  other  tanka, 
and  different  ambient  conditions.  Also,  Information  about 
the  transient  flow  dynamics  occurring  within  the  tank  was 
gleaned  from  the  heat  transfer  data.  This  Information 
Includes  rates  of  kinetic,  internal,  and  thermal  energy 
generation  and  dissipation  as  well  as  rates  of  entropy 
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UNSTEADY  HEAT  TRANSFER  RESULTING  FROM  THE  RAPID 

CHARGING  OF  AN  EVACUATED  TANK  WITH  CONDUCTING  WALLS 

I .  Introduction 

The  classical  problem  of  charging  a  tank  with  an  ideal 
gas  Is  well  documented  for  two  conditions  [1:80-871 
(2:203-2101.  If  an  Insulated  tank  Is  considered,  the  final 
equilibrium  temperature  can  be  determined.  If,  on  the  other 
hand,  the  tank  has  a  conducting  wall,  the  total  heat 
transferred  from  the  gas  to  the  wall  can  be  determined.  In 
either  case,,  the  solution  Is  derived  from  an  elementary 
application  of  the  first  law  of  thermodynamics.  The  change 
in  stored  energy  In  the  tank  must  equal  the  difference 
between  the  energy  which  came  Into  the  tank,  be  it  internal, 
kinetic,  or  flow  energy,  and  the  heat  lost  from  the  fluid  to 
the  tank  vails. 

In  contrast,  little  If  any  work  has  been  done  to 
understand  the  actual  transient  process  involved  In  charging 
a  vessel.  Thorough  literature  searches,  both  manual  and 
computer,  revealed  no  useful  Information  regarding  this 
process.  Perhaps  one  reason  for  the  lack  of  information  is 
that  no  analytical  solvitlon  exists  for  the  unsteady  process 
of  charging  a  tank  with  a  fluid.  Because  of  the  complexity 
of  the  flow,  an  analytical  solution  would  require  that  the 
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Navier-Stokes  equations  be  solved  in  their  entirety. 
Significant  simplifications  cannot  be  made  to  these  equations 
since  both  compreryulbi  1  i ty  and  diffusive  transport  of 
momentum  and  energy  dominate  the  entire  process.  The 
diffusive  processes  are  what  dictate  the  rate  of  kinetic 
energy  dissipation  in  the  tank  as  well  as  the  rate  of  heat 
transfer  to  the  walls.  As  a  consequence,  a  simple  energy 
balance  will  not  work.  The  only  alternative,  then,  is  an 
experimental  investigation.  With  the  rapid-response, 
thln-fllm  resistance  gages  currently  available,  the  rate  of 
heat  transfer  encountered  during  the  charging  of  a  tank  may 
be  obsorved  directly,  thereby  yielding  some  insight  into  the 
dynamics  of  the  flow. 

There  were  three  main  objectives  of  the  study.  The 
first  objective  was  to  detail  the  distribution  of  the 
measured  thermal  energy  along  the  wells  of  the  tank  by  direct 
observation.  The  second  abjective  was  to  use  the  measured 
heat  flux  data  to  derive  an  average  convection  coefficient 

-mJ  1  i.  n  I  m  ii  L  -  Aril  «  i>>  W  I  £  £  AMAMMJk  TU  AM  A  MMI  I  1  H  f  K  AS 

auu  i  V9  acouuviav,uu  tbiupcL  avvL  c  uiiicic.tuc.  auwoo  wu  «  %*  viivii 

be  used  to  predict  heat  flux  rates  encountered  during  the 
charging  of  other  tanks.  Finally,  the  third  objective  was  to 
gain  some  knowledge  of  the  flow  dynamics  that  occur  within  a 
tank  during  the  charging  process. 
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Charging  Is  the  process  by  which  a  vessel  containing  a 
fluid  at  low  pressure  is  filled  from  a  reservoir  of  fluid  at 
a  high  pressure.  The  resulting  pressure  gradient  forces 
fluid  from  the  reservoir  into  the  vessel  until  the  pressures 
inside  and  out  are  equal.  The  work  done  in  moving  this  mass 
of  fluid  is  then  manifested  in  the  internal  and  kinetic 
energies  of  the  fluid  in  the  vessel  as  well  as  any  heat  lost 
to  the  surroundings. 

For  example,  assume  that  a  perfectly  evacuated, 
insulated  tank  is  suddenly  opened  to  atmospheric  pressure. 
This  is  one  of  the  two  classical  problems  that  has  been 
solved.  The  mass  which  enters  the  tank  carries  with  it 
Internal  energy  as  well  as  flow  work  done  by  the  atmosphere 
in  pushing  the  mass  into  the  tank.  This  total  enthalpy  per 
unit  mass  is  Just  cpT0.  After  the  air  comes  to  rest  in  the 
tank,  it  has  energy  per  unit  mass  of  cvTf.  See  Figure  2-1. 

^  I  nr*o  fhn  fanb*  lex  rumrf  onf  1  \J  i  nail  1  AtnH  .  nri  psnf^mu  I  <1  1  nflt  tfl 


the  walls  by  way  of  heat  transfer.  Thus,  a  simple  energy 
balance  can  be  written  to  find  the  final  temperature  of  the 
air  In  tne  tank.  Enthalpy  in,  mcpT0,  must  equal  energy 
stored  in  the  tank,  mcyTj.  The  final  temperature  of  the  gas, 
therefore,  equals  T0cp/cv .  For  air,  if  TQ=*  530  R,  the  final 
temperature  would  be  742  R,  an  increase  of  212  R. 
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In  the  second  classical  problem,  the  tank  walls  are 
conductors  of  heat,  and  the  total  heat  transferred  to  these 
walls  is  sought.  Once  again,  assume  an  evacuated  tank  that 
is  opened  suddenly  to  the  atmosphere.  After  the  air  in  the 
tank  comes  to  equilibrium,  it  will  have  the  same  pressure  and 
temperature  as  it  had  before  entering  the  tank.  See  Figure 
2-2 .  All  the  flow  work  done  by  the  atmosphere  will  be 
manifested  in  the  heat  transfer  to  the  tank  walls.  This  flow 
work  per  unit  mass  is  Pq/^.  Therefore  Q  -  mPc/eo  -  PQV 
where  PQ  is  the  atmospheric  pressure  and  V  is  the  volume  of 
the  tank. 

During  the  transient  charging  period,  the  actual  rise  in 
fluid  temperature  within  the  uninsulated  tank  Is  a  complex 
function  of  the  rate  of  kinetic  energy  dissipation  in  the 
tank  as  well  as  the  rate  of  heat  dissipation  to  the  tank 
walls.  To  determine  the  heat  that  Is  lost  as  a  function  of 
time,  the  first  law  of  thermodynamics  in  its  open,  unsteady 
form  is  needed.  This  law  states  that  the  increase  in  stored 
energy  within  the  tank  must  equal  the  enthalpy  added  by  the 
mass  which  enters  the  tank  minus  the  heat  that  is  removed 
from  the  gas  at  each  instant  in  time.  Consider  the  control 
volume  pictured  In  Figure  2-3.  Mathematically, 
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Also, 

dE  »  d<  mcvT  +  mu^/2gc  )  [2-21 
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mcpT0  dt 


d<  rac„T  mu2/2ge  > 


Q  dt 


Figure  2-3  Control  Volume  Used  in  the  Transient  Energy 

Balance  of  Equation  (2-31 


»»o 


Figure  2-4  Control  Volume  Used  in  the  Transient  Entropy 

Balance  of  Equation  [2-41 


T  and  u  represent  some  average  temperature  and  velocity  of 
the  mass  within  the  tank. 

Equation  (2-23  can  be  substituted  Into  Eq.  [2*13  and, 
assuming  constant  enthalpy  entering  the  tank,  the  resulting 
Integrals  evaluated,  yielding 

m£cvT  +  u2/'2gc3  -  mjCvT0  -  <m-mj >cpT0  -  Q  C2-33 

The  term  rojcvT0  accounts  for  the  Initial  mass  In  the  tank  If 
It  was  not  fully  evacuated.  Q  is  the  total  heat  that  has 
been  transferred  from  the  gas  to  the  tank  walls. 

A  second  law  analysis  of  this  transient  problem  can  also 
be  accomplished.  The  second  law  of  thermodynamics  states 
that  the  rate  at  which  entropy  is  created  within  the  control 
volume  must  be  greater  than  or  equal  to  the  rate  at  which 
entropy  enters  with  the  entering  mass  minus  the  rate  at  which 
entropy  1-aves  by  way  of  heat  transfer.  See  Figure  2-4. 

This  can  be  expressed  as 

Sc,v.  -  "»o  ”  [2'4] 

Sc  v  Is  the  rate  of  entropy  Increase  within  the  tank,  sQ  Is 

e 

the  specific  entropy  of  the  fluid  entering  the  tank,  and  Q/T 
Is  the  rate  at  which  entropy  leaves  the  gas  due  to  heat 
transfer.  T  Is  the  temperature  of  the  Inner  tank  wall  which, 
for  walls  having  large  heat  capacity  and  conductivity,  equals 
Tq.  The  entropy  per  unit  mass  of  the  air  within  the  tank  can 
be  found  by 
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act)  -  aQ  +  cp  ln£T<t)/T0J  -  R  lnEP<t)/P03  £2-5] 

Thus,  data  derived  from  Eg.  £2-5}  can  be  numerically 
differentiated  to  flna  Sc  v  Ct>  providing,  of  course,  that 
the  mass  within  the  tank  is  known  since 

Sg ' y ' ( t )  -  d<  met)  act)  )/dt  £2-6] 
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HI  •  Thln-fllm  Gage  Theory 


Surface  Temperatures 

Thln-fllm  resistance  gages  were  used  to  measure  the 
temperature  at  the  surface  of  a  solid  substrate.  On  the 
surface  of  the  substrate  Is  a  thin  strip  of  platinum  through 
which  a  small,  constant  current  Is  passed.  The  resistance 
change  of  the  platinum  film  Is  directly  proportional  to  its 
change  In  temperature  above  some  reference  temperature.  The 
resistance  of  the  film  at  a  given  temperature  is 

Rg<T>  -  Rg<Tref>  +  <T-Tref>  <AR/AT>g  £3-13 

CAR/AT>g  Is  a  property  of  an  individual  gage  determined 
from  calibration.  C3»53 

The  resistance  change  of  the  film  can  be  measured  by 
measuring  the  voltage  change  across  it.  A  gage  can  be 
connected  to  a  bridge  and  amplifier  circuit  as  shown  in 
Figure  3.  The  bridge  element  resistor  that  is  shown  In 
series  with  the  gage  should  be  large  enough  to  ensure  that 
the  change  In  gage  resistance  does  not  induce  an  appreciable 
change  in  current.  For  this  type  of  constant-current  bridge 
circuit,  the  resistance  change  Is  given  by  the  measured 
voltage  change. 

A  Rg  ■  A  E  <R+Rg>/RIgG  £3-23 
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E  is  the  measured  voltage  across  the  bridge,  Ig  is  the 
constant  current  through  the  gage,  and  G  is  the  gain  of  the 
amplifier.  C31 

By  using  Eq.  1 3- 1 3  and  13-2],  the  change  in  voltage  can 
be  directly  converted  to  change  in  temperature  with  the 
following  expression. 

AT  -  A  E  <  R+Rg<Trof >+R'£T  -  Trftf]  >2  /  RVGR'  13-31 

where  R'  «  <AR/AT)g.  Note  that  current  has  been  eliminated 
in  favor  of  the  applied  voltage,  V. 


Surface  Heat  Fluxes 

Although  thin-fllm  gages  measure  temperature,  he&t  flux 
can  be  determined  by  using  a  solution  to  the  basic  heat 
equation  with  two  assumptions  Imposed.  First,  the  gage  is 
treated  as  a  semi-infinite  solid  during  the  test  time. 

Second,  the  heat  flux  through  the  gage  is  assumed  to 

prOpoQdtv  in  a  ufid”dl  nwHo  i  Ofid  1  fashion.  Tnuav  clcuicHiiS  a  i  i  OW 
heat  flux  rates  to  be  determined  from  the  temperature  data. 

To  solve  for  the  surface  heat  flux  in  terms  of  surface 
temperature,  the  one-dlmenslonal  heat  equation  can  be  used 

T+t  ■  a  T+xx  [3-41 

where  T+  *  T4<x,t>  -  T<x,t>  -  T<x,0>. 

The  boundary  conditions  are 
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T+CO,t>  -  T+s 
T+<x-»®,t>  -»■  0 
T+  <  x , 0  >  -  0 

The  solution  to  this  equation  can  be  found  by  using  a 
similarity  variable  or  by  the  the  Laplace  transform 
technique.  [4*384-390] 

T+  -  T+s  erfc<  x  /  £2<at>1/23  )  C3-5] 

An  energy  balance  at  the  surface  of  the  gage  requires  that 
the  heat  flux  incident  on  the  gage  must  equal  the  heat  flux 
into  the  gage.  Mathematical ly. 


q<t>  -  -  k  T+x!x-0 
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Duhamnel's  superposition  method  £4*403-4111  to  obtain  a 
numerical  solution  for  surface  heat  flux  given  a  discrete 
number  of  measured  surface  temperatures.  Their  equation, 
slightly  rearranged,  is 


q<tN> 


pck 

7T 


TV‘o> 

JZF* 


A _ +  _ ■ 

l  ”  1 


£3-73 


Thus,  given  the  measured  temperatures  from  a  thin-filra 
resistance  gage,  surface  heat  fluxes  may  be  calculated  as  a 
function  of  time. 
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IV.  The  Experiment 


Apparatus 

The  Tank.  The  tank  was  fabricated  from  a  section  of 
stainless  steel  pipe  0.375  Inch  thick.  The  tank  had  an  Inner 
diameter  of  10  Inches  and  a  length  of  35  13/16  Inches.  Holes 
were  bored  and  tapped  at  various  loc&tiono  on  the  tank  to 
accommodate  Instrumented  brass  plugs.  See  Figure  4-1.  The 
top  and  bottom  plates  were  made  from  steel  and  provided  with 
sites  for  instrumented  plugs  as  well.  The  bottom  plate  was 
welded  in  place,  while  the  top  one  was  bolted  to  the  tank. 
This  provided  a  means  of  entering  the  tank  should  the  need 
arise.  All  touching  metal  parts  that  were  not  welded 
together  had  G-rlng  seals  or  rubber  gaskets  between  them  to 
prevent  leakage. 

The  entrance  to  the  tank  consisted  of  a  simple 
converging  nozzle.  The  nozzle  had  an  entrance  diameter  of  4 
inches  and  a  throat  diameter  of  1  inch,  yielding  an  area 
ratio  of  16.  A  rubber  stopper  commonly  found  in  chemistry 
labs  was  used  to  plug  the  nozzle  during  evacuation.  The 
stopper  was  attatched  to  one  arm  of  a  lever  which  was  bolted 
to  the  top  of  the  tank.  This  lever  was  used  to  initiate  the 
flow  quickly  by  pulling  the  stopper  out  and  away  from  the 
nozzle  Inlet. 

Th in- film  Gages .  Five  thln-fllm  gages  were  used  to 
gather  data  on  the  rates  of  heat  transfer  at  19  different 
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stations  over  the  taak.  The  gages  were  commercially  produced 
and  consisted  of  a  thin  strip  of  platinum  which  was  deposited 
onto  one  end  of  a  cylindrical  piece  of  quartz.  The  quartz 
substrate  was  0.25  inch  long  and  had  a  diameter  of  0.125 
Inch.  Each  gage  was  Inserted  Into  a  hollow  Plexiglas  rod 
which  had  an  outer  diameter  of  0.3875  inch.  This  rod  was 
then  affixed  inside  a  hollow  brass  plug.  See  Figure  4-2. 
Epoxy  was  used  to  hold  the  plexlglas  and  quartz  in  place  as 
well  as  to  ensure  an  airtight  seal.  When  Inserted  into  the 
tank,  the  center  of  each  gage  was  flush  with  the  inner  wall 
of  the  tank;  the  edges  were  slightly  recessed  due  to  the 
curvature  of  the  tank  wall. 

The  gages  were  connected  to  br 1 dge/amp 1 1 f 1 er  circuits. 
The  output  from  each  bridge  was  stored  using  a  digital 
waveform  recorder.  After  a  run,  each  of  the  five  recorded 
channels  was  displayed  on  an  osc i 1 1 1  scope ,  and  a  Polaroid 
photograph  was  taken  of  the  resulting  trace.  The  photo  was 
later  digitized  and  the  data  stored  in  a  computer  for  data 
reduction. 

Pressure  Transducer .  An  Endevco  0-100  psia  pressure 
transducer  was  used  to  measure  the  transient  pressure  within 
the  tank.  The  transducer  and  its  signal  conditioner  fed  data 
to  the  same  digital  waveform  recorder.  A  Polaroid  photograph 
was  taken  and  digitized  and  the  data  was  stored  in  a  computer 
for  later  reduction.  A  mercury  manometer  was  used  to 
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calibrate  the  transducer  while  a  mercury  barometer  was  used 
to  determine  atmospheric  pressure. 

Vacuum  Pumo .  A  mechanical  vacuum  pump  was  used  to 
evacuate  the  tank.  The  pump  was  connected  to  a  valve  which 
was  attached  to  the  side  of  the  tank.  It  was  consistently 
able  to  reduced  the  tank  pressure  to  between  2  and  3  psia 
within  5  to  7  minutes. 

Digital  Waveform  Recorder.  A  Datalab  8  channel  digital 
waveform  recorder  was  used  to  capture  the  data  from  the 
thin- film  and  pressure  gages.  It  was  desired  to  use  the 
digitized  data  directly,  but  a  suitable  computer  was 
unavailable.  Thus,  the  information  stored  in  the  recorder 
was  displayed  on  an  osc 1 1 1 1  scope ,  photographed,  and  then 
digitized  using  other  equipment. 

Dial t izer/Comouter ♦  A  Hewlet-Packard  plotter/  digitizer 
was  used  to  digitize  the  thin-film  and  pressure  traces.  This 
data  was  then  input  by  hand  into  an  Apple  lie  computer  for 
reduction  and  manipul at  Ion . 
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Pressure  Gage  Calibration 

To  calibrate  the  pressure  transducer,  the  tank  was 
evacuated  eight  different  times  to  various  pressure  levels. 
Each  time,  the  pressure  from  a  mercury  manometer  was  recorded 
along  with  the  deflection  of  the  transducer  trace  on  the 
ocsl 1 1 1  scope ,  These  measurements  gave  conversion  constants 
(from  millivolts  to  psl)  which  varied  by  not  more  than  4%. 
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These  eight  constants  were  then  averaged  together  to  give  a 
conversion  constant  of  2.88  mV/psl . 

Next,  to  determine  which  location  on  the  tank  would 
give  the  "best*  average  pressure,  the  transient  pressure  was 
recorded  at  the  top,  the  middle,  and  the  bottom  of  the  tank. 
There  was  no  perceptible  difference  in  the  three  traces.  As 
a  result,  the  pressure  was  assumed  to  be  essentially  constant 
throughout  the  entire  tank  at  each  Instant  in  time.  Since 
the  location  of  the  transducer  was  unimportant,  the  middle 
location  was  chosen  for  all  subsequent  runs  as  a  matter  of 
convenience. 


Heat  Flux  Gage  Calibration 

Due  to  the  difficulty  in  calibrating  these  types  of 
gages  (33,  the  calibration  from  the  manufacturer  was  assumed 
to  be  correct.  Each  of  the  gages  had  different  resistances 
and  different  rates  of  resistance  change  with  temperature 
change.  Care  was  taken  during  all  test  runs  to  keep  track  of 
each  gage  and  the  position  on  the  tank  which  it  occupied  at 
that  time. 

Because  thln-fllm  gages  work  on  a  one-dlmensionai  heat 
flow  model  (see  Section  III),  the  long  times  of  this 
experiment  (l.e.  >  1  second)  caused  some  concern  as  to 
whether  two-dimensional  effects  might  al -er  the  results.  In 
order  to  find  out,  one  of  the  gages  was  exposed  to  a  constant 
heat  flux  and  the  resulting  signal  recorded.  Since  an 
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analytic  solution  exists  for  this  boundary  condition,  the 
response  of  the  gage  could  be  compared  to  a  true  one- 
dimensional  model.  The  source  of  the  constant  heat  flux  was 
the  stagnation  region  of  a  hot  stream  of  air.  This  stream 
was  generated  with  a  hot  air  gun  which  is  normally  used  to 
shrink  electrical  insulation.  In  actuality,  it  is  the 
Mussel t  number  or  convection  coefficient  which  remains 
constant  in  stagnation  region  heat  transfer.  However,  if  the 
surface  temperature  of  tho  gage  is  assumed  to  remain 
essentially  constant,  then  heat  flux  will  also  be  constant. 

The  test  was  crude  and  the  hot  stream  only  somewhat 
constant,  but  the  results  Indicated  that  the  gage  was  giving 
approximately  one-dimensional  results  all  the  way  out  to  ? 
seconds.  According  to  the  analytic  solution,  a  constant  heat 
flux  would  result  in  an  increase  in  surface  temperature 
proportional  to  the  square  root  of  time.  This  would  generate 
a  horizontal  line  as  plotted  in  Figure  4-3.  Between  0  and  1 
second  for  the  low  heat  flux  levels,  there  were  small 
deviations  from  a  horizontal  line.  This  was  due  to  the 
transient  effects  associated  with  establishing  the  stream  of 
air  onto  the  gage.  The  difference  was  much  greater  for  the 
two  highest  heat  flux  levels  where  It  was  necessary  to  have 
the  fan  on  the  gun  set  on  high.  Apparently,  this  Increased 
the  transient  flow  establishment  time.  After  the  transient 
time  period,  however,  the  curves  do  approximate  a  straight, 
horizontal  line.  However,  the  curves  generated  with  the  two 


Figure  4-4  Heat  Flux  Curves  Used  in  Determining 

Exper iroent  Repeatabi l 1 ty 


heat  flux  levels  seem  to  continue  decreasing  after 
the  Initial  transients.  Because  these  levels  produced  an 
appreciable  surface  temperature  change  (near  10  F) ,  the  heat 
flux  was  no  longer  a  constant  and  was  expected  to  decrease. 

As  a  result,  no  correction  was  employed  to  account  for 
two-dimensional  effects. 

To  ensure  that  the  main  experiment  was  repeatable,  three 
preliminary  runs  were  made  with  the  gages  mounted  in  the 
tank.  The  results  of  one  of  the  gages  is  shown  in  Figure  4-4. 
The  three  traces  were  nearly  Identical:  the  slight  variations 
due  largely  to  the  slightly  differing  initial  pressures  in 
the  tank.  Although  convinced  of  the  repeatability  of  the 


\  tiuvin*  | 


Cle  w  W«S  Se«  1  18  VW  A  HVW 


1 l aa 
■  vvw  w  4  wii 


throughout  all  subsequent  runs  to  ensure  that  this  assumption 
remained  valid. 
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before  a  test  was  n?ade;  the  electronic*  were  alven  at 
least  20  minutes  to  warm  up,  allowing  all  transients  to  die 
out.  The  tank  was  then  evacuated  to  the  capacity  of  the 
vacuum  pump:  this  was  between  2  and  3  rsia.  During  the 
evacuation,  atmospheric  pressure  and  temperature  were 
recorded.  After  that,  the  bridge  circuits  were  balanced,  and 
the  digital  recorder  was  checked  to  ensure  that  the  settings 
employed  were  the  correct  ones.  The  first  five  channels  of 
the  recorder  were  used  for  the  five  heat  flux  gages,  the 
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sixth  for  the  pressure  transducer,  and  the  remaining  two  went 
unused.  The  recorder  was  set  to  trigger  from  any  change  in 
the  pressure  tranducer  since,  from  earlier  testa,  it  seemed 
to  have  to  fastest  response.  The  digital  recorder  was  also 
set  to  record  data  for  6  seconds. 

Once  the  tank  had  been  evacuated,  the  valve  between  it 
and  the  pump  was  closed;  then,  the  pump  was  turned  off.  The 
recorder  was  armed  and  the  tank  stopper  lever  was  pulled  down 
sharply  to  ensure  a  rapid  openning.  After  the  8  seconds  were 
over,  each  of  the  recorded  channels  was  displayed  on  an 
osci 1 1 1 scope.  The  image  was  then  enlarged  and  expanded  using 
the  controls  on  the  recorder.  This  was  done  to  facilitate 
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tests  Indicated  that  the  significant  heat  transfer  occurred 
within  the  first  4  seconds,  only  this  portion  of  the  signal 
was  photographed.  To  take  the  pictures,  the  grid  on  the 
oscilliscope  was  illuminated  and  the  intensity  of  the  trace 
made  slightly  brighter  than  the  grid.  The  camera  was  set 
with  a  shutter  speed  of  1/30  second  and  an  f-stop  of  4.  Each 
photograph  was  marked  with  an  identifying  code  immediately 
after  processing.  This  code  was  used  to  identify  the  gage 
and  tank  position  for  that  particular  trace.  Figures  4-5  and 
4-6  are  sample  photographs  of  a  heat  flux  trace  and  pressure 
trace,  respectively. 
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At  a  later  date,  the  photographs  were  digitized  using  a 
Hewlet-Packru  d  plotter/digitizer.  Because  the  computer  to  be 
used  in  the  daca  reduction  could  not  be  made  to  communicate 
with  the  digitizer,  data  from  the  digitizer  had  to  be  entered 
manually  into  the  computer.  This  was  the  most  laborious  task 
of  the  experiment  and  was  also  a  driving  factor  in  limiting 
the  recorded  signal  time  to  4  seconds. 
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Figure  4-5  Sample  Oac i 1 1 1 scope  Output  from  Th ln-fllm  Gage 

$ 


Figure  4-6  Sample  Oscllllscope  Output  from  Pressure 

Transducer 
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V.  Data  Reduction 


general  Caaagala 

From  the  five  different  tests  needed  to  accumulate  data 
from  all  19  locations  on  the  tank,  average  initial  conditions 
were  calculated  and  used  In  all  data  reduction. 

The  initial  conditions  used  were 

P0  »  2066.65  lbf/ft2  0.6  % 

P.  -  200.365  lbf/ft2  1.3  % 

T0  «  540  R  0.9  * 

The  numbers  in  the  far  right  column  Indicated  the  largest 

percentage  difference  from  the  average. 

In  all  subsequent  calculations,  numerical  derivatives 

were  calculated  using  a  centered  difference  formula  for 

interior  points  or  a  three-point  formula  for  end  points.  For 

any  integrations,  cubic  splines  were  fit  to  the  data  points 

and  the  resulting  cubic  polynomials  integrated  exactly. 

Measured.  Pat.a 


One  of  the  first  quantities  calculated  from  the  measured 
data  was  the  mass  in  the  tank  as  a  function  of  time.  The 
only  data  needed  for  this  calculation  was  the  transient 
pressure  data  and  the  initial  conditions.  Because  the 
entrance  to  the  tank  was  a  simple  converging  nozzle, 
lsentropic  relations  could  be  used  to  calculate  the  mass  flow 
as  a  function  of  the  pressure  In  the  tank.  In  terms  of  Mach 
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number  at  the  throat,  the  mass  flow  through  the  nozzle  is 
given  by  16:84-85] 

TO  *  PQ  A  M  <r/RT0>1/2  £  1  +  M2<Y-l>/2  3<v+1  )/2<  i~y>  £5-13 

For  a  perfect  gas  under  lsentropic  conditions,  the  Mach 
number  can  be  replaced  with  the  following  expressions. 


For  unchoked  flow, 

M  **  £2/<Y-l>31/2  £(P/P0><1-7>/'y-131/2 

For  choked  flow, 

M  -  1 


£5-23 


P  is  the  pressure  within  the  tank  and  Pn  is  the  atmospheric 
pressure.  Substituting  £5-23  into  £5-13  and  letting  P+=* 
CP/P0)(1~7)/y  yields  the  following. 


For  unchoked  flow, 

m  «  PQ  A  £2  /av-13RT0)31/2  £P+-13l/2  P+(7+1  >/2Cl-y> 


For  choked  flow, 

m  »  PQ  A  <Y/RT0)1/2  £  <7+13/2  jO'+l  >/2<  1-T> 


£5-33 


According  to  the  measured  pressure  data,  see  Figure  5-1,  the 
nozzle  unchoked  after  0.18  second,  and  the  pressure  in  the 
tank  reached  atmospheric  pressure  after  0.46  second. 

Equation  £5-33  was  Integrated  numerically  from  t*0  to  t=0.46 
second  to  obtain  the  mass  in  the  tank  as  a  function  of  time. 
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The  result  is  shown  In  Figure  5-2.  It  should  be  noted  that 
the  mass  In  the  tank  at  0.46  second,  0.11017  lbm.  Is  somewhat 
less  than  the  mass  the  tank  could  hold  at  equilibrium, 

0.11749  lbm.  This  Indicated  that  the  temperature  of  the  gas 
was  still  above  Tc . 

From  the  pressure  data  and  the  newly  calculated  mass 
data,  the  average  temperature  of  the  air  within  the  tank  was 
found  as  a  function  of  time.  This  temperature  was  found  by 
using  the  ideal  gas  law 

T< t >  »  P<t>  V  /  m<t>  R  15-41 

The  results  are  plotted  in  Figure  5-3.  Because  the  partial 
derivative  of  temperature  with  respect  to  mass  is 
proportional  to  the  reciprocal  of  the  mass  squared,  small 
errors  in  calculating  the  mass  were  magnified  when  the 
temperature  was  calculated.  This  explains  the  somewhat 
erratic  nature  of  Fig.  5-3.  However,  the  general  trend  of 
increasing  temperature  is  readily  apparent. 

Next,  from  the  voltage  data  of  the  thln-film  gages, 
surface  temperatures  were  calculated  using  Eq.  [3-31.  These 
temperatures  were  then  used  in  the  Cook  and  Felderman 
equation  [3-71  in  order  to  determine  the  surface  heat  fluxes. 
The  value  for  pck/'7r  used  in  this  equation  was  the  well 
documented  value  of  0.001819  BTU^/f t4-R-s.  13]  This  data 
spanned  the  time  period  between  0  and  4  seconds.  The  heat 
flux  of  each  gage  was  calculated  separately,  and  then  It  was 
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Figure  5-2  Maaa  of  Air  Within  Tank:  0-0.46  Second 


Figure  5-3  Average  Air  Temperature  in  Tank:  0-0.46  Second 
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mu  1 t 1 p 1 yed  by  the  area  of  that  section  of  tank  it  was 
representing.  The  results  from  all  the  gages  were  then  added 
together  to  get  the  total  heating  rate  for  the  entire  tank. 
See  Appendix  B  for  this  information  in  tabular  form.  The 
aata  was  then  integrated  numerically  to  determine  the  heat 
transferred  to  the  tank  walls  as  a  function  of  time.  This 
data  is  displayed  in  Figure  5-4. 

At  the  end  of  the  4  seconds,  che  total  heat  transferred 
was  approximately  80%  of  the  total  amount  which  should  have 
been  transferred  had  equilibrium  been  reached.  That  maximum 
is  3.92  BTU  and  depends  solely  on  the  initial  and  final  tank 
pressures  as  well  as  the  volume  of  the  tank.  (See  Section 
II:  the  second  classical  piobieiii.)  A»  a  result,  the  Curve 
was  extrapolated  out  to  10  seconds  to  account  for  the  missing 
energy,  but  that  will  be  discussed  shortly.  The  heat  flux 
data  from  the  individual  gages  versus  their  position  in  the 
tank  are  plotted  in  Figure  5-5  at  specific  moments  in  time. 
These  plots  give  an  indication  of  the  heating  distribution 
along  the  tank  wall. 

Fran  the  mass,  temperature,  and  heat  data,  the  average 
Internal  and  kinetic  energies  of  the  gas  within  the  tank  were 
calculated.  Recall  Eq.  [2-3], 

mCcvT  +  u^/2gc3  -  mjCvT0  =  (m-mj )cpT0  -  Q 
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Heat  Flux  Versus  Wall  Location 


Since  the  only  unknown  in  Eq.  12-3)  Is  the  kinetic  energy 
term,  mu2/2gc>  it  could  be  determined.  Each  of  the  energy 
terms  appearing  in  Eq.  C2-33  is  plotted  in  Figure  5-6. 

Because  rapid  charging  Involves  high-speed  flow,  the 
adiabatic  wall  minus  wall  temperature  difference  was  chosen 
as  the  driving  temperature  potential  for  heat  transfer.  Due 
to  the  unknown  dynamics  of  the  flow,  a  recovery  factor  of 
unity  was  chosen  for  lack  of  a  better  estimate.  This  reduced 
the  driving  potential  to  be  the  total  temperature  minus  the 
wall  temperature.  The  total  temperature,  T* ,  is  shown  in 
Figure  5-7. 

T*<t>  -  T< t )  +  tu(t>)2/2gccp  [5-53 

To  determine  an  average  convection  coefficient  for  the 
tank,  h(t),  the  following  equation  was  used. 

h<  t )  -  q<  t  >  /  LT*( t )  -  T03  [5-63 

To  get  qCt),  the  data  for  GKt)  was  divided  by  the  total 
surface  area  of  the  tank  and  numerically  dl f ferentlatied. 

The  plot  of  h<t)  is  shown  in  Figure  5-8.  Anticipating  a 
relationship  between  h<t>  and  Reynolds  number,  logCh)  was 
plotted  versus  the  log<  pu).  The  velocity,  u,  was  determined 
from  the  kinetic  energy  data.  This  logarithmic  plot  Is  shown 
in  Figure  5-9. 
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5-6  Terms  In  the  Energy  Balance  of  Eq.  [2-33* 

0-0.46  Second 


Ealcaefll  ated-Data 


Because  the  heat  transferred  after  4  seconds  was  short 
of  the  theoretical  equilibrium  value,  an  attempt  was  made  to 
account  for  this  energy  by  extending  the  time  for  equilibrium 
to  be  established.  After  several  attempts,  10  seconds  was 
found  to  be  the  time  at  which  the  heat  transfer  was  completed 
and  at  the  maximum  value.  In  order  to  extrapolate  all  the 
data  out  to  10  seconds,  three  quantities  had  to  be  chosen  as 
Independent  variables.  Since  the  pressure  within  the  tank 
wao  known.  It  was  selected  as  one  of  the  independent 
variables.  The  two  other  quantities  finally  chosen  were  the 
convection  coefficient,  h,  and  the  velocity,  u. 

The  convection  coefficient  was  chosen  as  an  Independent 
variable  because  prior  trials  at  extrapolation  indicated  that 
h<t>  remained  fairly  constant  after  1  second.  Therefore,  it 
was  assumed  that  from  1  to  10  seconds  h(t)  *  0.0014386 
BTU/ft^-F-s  .  This  modified  curve  Is  displayed  in  Figure 

c-  «  n 

The  velocity  was  also  chosen  as  a  known  quantity  because 
there  was  sow©  theory  with  which  tc  predict  its  behavior.  The 
rate  of  momentum  dissipation  of  a  fluid  within  a  tank  depends 
upon  the  force  or.  that  fluid  due  to  the  shear  stress  at  the 
wall.  Since  the  mass  changed  very  slowly  after  0.46  second, 
it  may  be  approximated  as  a  constant.  Doing  so,  the  momentum 
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age  Convection  Coefficient 
-10  Seconds 


dissipation  may  be  expressed  as 


d<  mu<t>  >/dt 


"  tw 


A 


surface 


After  integration, 


-  ^ u< t )Agur^ace/L  15-71 


u<t)  «  u(0.46)  exp  <-(  H  Agurjace/mL)(t-0.46)>  C5-83 

Where  rw,  wall  shear  stress,  has  been  approximated  by 
fxu(t.)/L,  and  L  is  some  characteristic  length.  From  Eq. 
[5-81,  the  velocity  should  decrease  exponentially.  Because 
the  velocity  prior  to  0.46  second  was  calculated  from 
measured  quantities,  a  simple  exponential  decay  was  fit  to 
the  velocity  curve  from  0.46  to  10  seconds.  The  time 
constant,  i/<  n  Aaurface/m£<> »  wa8  taken  to  b®  0.5  ®*c  in  order 
to  have  u(10)  approximately  equal  to  zero.  This  curve  is 
shown  in  Figure  5-11. 

To  find  the  mass  In  the  tank  as  a  function  of  time,  the 
energy  equation,  Eq.  C2-3],  can  be  rewritten  In  terms  of  P, 
h,  m,  and  u. 


cvPV/R  +  mu2/2gc  -  mjCvT0  * 

mcpT0  -  hA<PV/mR  +  u2/2gccp  -  T0>  15-9] 

Differentiating  with  respect  to  time  and  solving  for  m  yields 

<  muu/gc  )  +  dthA<PV/mR  +  u2/2gccp  -  TQ>]/dt 

m  =»  - - - - - -  15-10] 

CpTo  -  u2/2gc 

Equation  [5-10]  was  numerically  integrated  and  the  resulting 
mass  function  is  plotted  in  Figure  5-12.  With  the  mass  in 


the  tank  now  estimated,  the  average  and  total  temperatures 
were  found  using  Eqs.  £5-4]  and  £5-5],  respectively.  They 
are  shown  in  Figure  5-13.  From  this  resulting  data,  the 
energy  terms  in  Eq.  £2-3]  were  calculated.  They  are 
displayed  in  Figure  5-14.  Also,  the  plot  of  log(h)  versus 
log(pu)  was  found  and  is  shown  in  Figure  5-15. 

To  verify  that  the  second  law  of  thermodynamics  was  not 
violated  as  a  result  of  the  modified  data,  each  of  the  three 
terms  of  Eq.  £2-4]  was  calculated  from  the  estimated  data. 

As  a  convenience,  entropy  per  unit  mass  was  taken  to  be  zero 
at  atmospheric  conditions.  This  eliminated  the  need  to 
calculate  the  ms0  term. 

The  tank  wall  temperature  was  taken  to  be  T0  since  the 
wall  had  such  a  large  heat  capacity  and  the  heat  flux  was  so 
small.  In  fact,  for  a  constant  heat  flux  applied  to  the 
boundary  of  a  semi-infinite  solid,  the  ratio  of  temperature 
change  for  a  stainless  steel  solid  (tank  wall)  to  that  of  a 
quartz  solid  <th ln-fllm  gage)  is  given  by 

i^Ttank/  ^Tgage  “  c  Pck)gage/<  Pck)tank  £5-113 

For  the  material  properties  involved,  this  ratio  is 
approximately  0.2.  Since  the  maximum  gage  temperature 
differences  were  between  3  and  5  F,  the  maximum  temperature 
change  of  the  tank  wall  would  be  about  1  F.  Thus,  the  wall 
temperature  was  assumed  to  remain  constant  at  TQ. 

To  calculate  Sc  v  ,  Eq.  £2-5]  was  used. 
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Extrapolated  Average  and  Total  Temperatures 
0-10  Seconds 


p  1 n t T/T0 3  -  R  lnEP/P0l 


[5-123 


s<t> 


Thus, 


sc.v,(t)  *  m<t)s<t> 


or. 


[5-133 


Sc  v  <t>  *  d<  ra< t )s< t)  >/dt 


The  results  of  these  calculations  are  shown  in  Figure  5-16. 
Recall  that  sc  *  C.  Equation  [2—43  then  requires  that 


sc.v.<t>  ♦  q/td  l  o 


[5-143 


Figure  5-1?  shows  the  sum  of  these  entropy  terms,  and  indeed, 
the  sum  is  initially  greater  than  and  quickly  approaches 


zero. 


Entropy  Rates  <  BTU/s- 


Entropy  (  BTU/F 
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VI. 


Results 


As  shown  In  Figure  5-5,  the  heat  flux  level  along  the 
tank  wall  at  any  given  time  is  fairly  uniform.  Deviations  do 
exist,  however,  and  may  Imply  recirculation  regions  within 
the  tank.  With  the  exception  of  the  positions  on  the  bottom 
of  the  tank,  the  crests  In  heat  flux  remain  crests  and  the 
troughs  remain  troughs.  What  may  be  happenning  in  the  tank 
is  this:  the  Instant  the  tank  is  opened,  air  flowing  out  of 
the  choked  nozzle  expands  Immediately  into  the  upper  corners 
of  the  i.ank.  The  flow  then  continues  down  the  tank  and 
finally  stagnates  at  the  base.  As  the  tank  fills,  the  flow 
out  of  the  nozzle  becomes  more  jet-like  and  swirling  begins. 
This  swirling  seems  to  reach  a  peak  around  the  time  at  which 
the  nozzle  unchokes.  After  that  time,  the  velocity  of  the 
entering  air  decreases  as  does  the  kinetic  energy  of  the  air 
within  the  tank.  The  flow  pattern  during  the  time  of  maximum 
swirling  may  appear  as  that  shown  in  Figure  6-1. 

From  all  the  figures  in  Section  V,  it's  clear  that  the 
bulk  of  the  charging  process  is  over  in  a  very  short  time. 

The  pressure  within  the  tank  rose  to  atmospheric  pressure  in 
0.46  second.  The  mass  nearly  rose  to  its  maximum  value  in 
0.46  second.  From  that  time  on,  mass  slowly  entered  the  tank 
as  the  average  gas  temperature  within  gradually  decreased. 
After  0.46  second,  slightly  over  50%  of  the  heat  energy  had 
already  been  transferred  from  the  gas  to  the  tank  walls.  By 
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four  seconds,  80%  had  been  transferred;  by  8  seconds,  95%; 
and  by  10  seconds,  almost  100%. 


A  convection  coefficient  was  derived  and  plotted  in 
Figure  5-10.  This  h(t)  was  also  plotted  on  a  log  scale  and 
Is  shown  in  Figure  5-15.  The  plot  contains  the  same  strange 
behavior  as  the  plot  In  Fig.  5-9.  This  data  Indicates  a 
complex  flow  pattern  within  the  tank  with  two,  possibly 
three,  different  regimes.  During  the  early  times  of  the 
process,  h  Is  proportional  to  (pu)1-***.  A  transition  then 
occurs  and  h  seems  to  be  independent  of  <pu>.  The  third  flow 
pattern  is  the  straight  line  in  the  plot  and  is  a  result  of 
choosing  h(t>  to  be  constant  after  1  second.  Little  can  be 
said  concerning  the  meaning  of  that  portion  of  the  curve. 
Because  of  this  strange  behavior.  It  waa  not  possible  to 
establish  a  relationship  between  h<t>  and  Reynolds  number. 

This  erratic  behavior  of  the  convection  coefficient  also 

precluded  its  use  as  a  tool  for  predicting  heating  rates  in 

other  tanks.  Another  approach  was  needed.  To  this  end,  the 

« 

neat 1 1 1 y  rate  for  the  entire  tank.  Gets,  was  found  from 
differentiating  GKt)  and  is  displayed  in  Figure  6-2. 

0<t)  depends  upon  the  initial  pressure  conditions,  the 
volume  of  the  tank,  and  the  size  of  the  inlet.  Also,  it  will 
depend  upon  the  time  during  which  the  nozzle  is  choked  and 
the  time  for  pressure  within  the  tank  to  equalize  with  the 


ambient  pressure.  These  factors  were  considered  in  extending 
the  measured  results  to  other  tanks  with  conducting  walls 


which  have  roughly  the  same  thermal  properties  as  stainless 
steel.  Those  properties  are  a  =  0.147  ft2/hr  and  pck  = 

435  BTU2/f t4-hr-F2.  It  should  be  noted  here  that  the 
following  calculations  (with  the  exception  of  Qtot>  are  only 
estimates. 

As  was  discussed  previously,  the  total  heat  transferred 
from  the  gas  to  the  tank  walls  is  just  the  total  flow  work 
done  on  the  gas.  This  energy,  Qtot,  equals  the  difference  in 
initial  pressures  multiplied  by  the  volume  of  the  tank.  This 
relationship  is  plotted  in  Figure  6-3  for  various  heat 
levels. 

To  find  the  time  at  which  the  nozzle  unchokes  as  well  as 
the  time  at  which  the  pressures  equalize,  the  measured 
pressure  curve  was  represented  by  two  straight  lines.  The 
first  line  describes  the  pressure  rise  while  the  nozzle  is 
choked.  Since  this  relationship  is  a  linear  one,  no 
approximation  needs  to  be  made.  The  second  line  represents 
the  pressure  rise  between  the  time  from  the  unchoklng  of  the 
nozzle  to  the  time  the  tank  pressure  reaches  PQ.  This  is  an 
approximation  to  the  actual  curve.  See  Figure  6-4.  The 
slopes  of  these  lines  were  calculated  from  the  measured  data, 
and  they  yielded  the  following  results. 

tuc  -  <0.52823Po  -  Pp  /  4952.25 

[6  -13 

tp  =  tuc  +  U  -  0.52828)  PQ  /  3481.71 
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6-3  Relationship  Between  Pressure  Diffe 
Tank  Volume,  and  Total  Heat  Transf 


■‘M*,  The  pressures  are  measured  in  lbf/ft^,  and  the  times  are  in 

seconds.  These  two  equations  are  only  valid  for  the  nozzle 
and  tank  that  was  used  in  this  study.  Since  the  time  to  fill 
a  tank  ought  to  be  proportional  to  the  volume  of  the  tank, 

Eq.  [6-1]  was  multiplied  by  the  ratio  of  new  tank  volume,  V, 
to  the  tank  volume  used  In  this  study.  Also,  a  nozzle  with  a 
large  throat  area  will  fill  a  tank  faster  than  one  with  a 
small  throat  area.  Thus,  Eq.  £6-1]  was  divided  by  the  ratio 
of  new  throat  area,  Ath,  to  the  area  of  the  throat  that  was 
used.  These  equations  are 


tuc  «  0.003362  <V/AthX0.52828Po  -  Pi  >  /  4952.25 

[6-23 

tp  =  tuc  +  0.001585  <VP0,'Ath>  /  3481.71 

The  tank  volume  is  measured  in  ft3  and  the  throat  area  in 
ft*.  Equation  [6-23  can  now  be  used  to  estimate  the  time  it 
takes  for  a  9iven  nozzle  to  unchoke  as  well  as  the  time 
required  for  the  gas  in  the  tank  to  reach  equilibrium 
pressure. 

Frcsm  the  data  in  Figure  5-14,  the  the  heat  transferred 
by  tp  is  approximately  50%  of  0to£.  Also,  from  Figure  6-2, 
CKt>  can  be  accurately  approximated  by  a  second  order 
polynomial  between  t*0  and  t»tp. 

QCt>  =  [at2+  bt  +  c 3  [6-33 
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The  following  equation  can  then  be  used  to  find 

0.5  Qtot  -  ^x  J  tat2+  bt  +  cl  dt  [6-41 

The  three  conditions  needed  to  find  tho  coefficients  are 

Qmax  *  ^Ctucj 
QCO)  =  0 

6<tp>  -  0.18  Qkm 

Therefore  , 

OCt)  *  Gfoaxt  t8.875  -  18.444  tl  [6-5] 

After  substituting  into  Eq.  C6-41  and  Integrating, 

Qjjjax  ®  0.5  Qtot  /  t  tp2<4 . 4375  -  6.i48tp)3  (6-6i 

Thus,  by  knowing  only  the  volume  of  a  tank,  the  area  of 
the  nozzle  throat,  the  ambient  pressure,  and  the  initial 
pressure  within  the  tank,  the  total  heat  transferred  to  the 
tank  walls  can  be  calculated.  Also,  the  times  to  unchoklng 
and  pressure  equalization  can  be  estimated  as  can  the  peak 
heat  flux  and  the  shape  of  the  GCt)  curve.  The  remainder  of 
the  QCt)  curve  can  also  be  estimated  after  tp  by  assuming 
another  second  order  polynomial .  This  part  of  the  curve 
represents  long  times  with  small  heating  rates.  Because  of 
this,  the  first  part  of  the  curve  is  most  important  and  was 
all  that  was  found.  However,  the  total  time  for  charging  to 
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be  completed  can  be  estimated  by 


tf  =  22  tp  C6-6] 

This  equation  was  found  by  taking  the  measured  ratio  of  the 
total  time  for  charging,  10  seconds,  to  the  time  at  which  the 
tank  pressure  equaled  PQ,  0.46  second. 
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vii  .  summac*. 


This  study  was  a  new  look  Into  the  old  problem  of 
charging  a  vessel  with  an  ideal  gas.  The  actual  transient 
process  was  observed  and  calculated.  A  major  suprise  was  the 
uniformity  of  heat  flux  along  the  walls  of  the  tank.  There 
were  no  significantly  hot  or  cold  spots  at  any  point  on  the 
tank  wall.  Another  suprise  was  the  strange  behavior  of  the 
convection  coefficient  with  respect  to  Cpu).  As  a  result,  it 
was  difficult  to  make  any  conclusive  statements  as  to  the 
actual  flow  dynamics.  It  was  also  difficult  to  formulate  a 
method  of  using  h<t>  for  the  prediction  of  heat  transfer 
rates  in  other  tanks.  However,  a  method  of  estimating  peak 
heating  rates  as  well  as  the  general  shape  of  the  Q<t)  data 
was  found. 
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VIII.  Recommendations 


Further  study  of  the  charging  process  is  definitely 
warranted.  This  study,  being  the  first  of  it's  kind,  has 
only  scratched  the  surface  of  the  complexities  Involved  in 
the  flow  dynamics  and  heat  transfer  during  the  charging  of  an 
evacuated  tank.  Further  experimentation  should  involve  the 
recording  of  heat  transfer  data  throughout  the  entire 
process.  Also,  testing  with  different  nozzle  sizes  and 
various  pressure  conditions  should  be  carried  out  to  verify 
or  disprove  the  est inflation  techniques  developed  here.  It 
would  also  be  of  Interest  to  vary  the  length  to  diameter 
ratio  of  the  tank  in  order  to  determine  effects  resulting 
from  geometry.  All  these  would  aid  In  developing  a  means  of 
accurately  predicting  heat  flux  during  the  charging  process. 
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Table  of  Heats  Measured  by  Individual 
Thln-Fl 1m  Gages 


*  LOCATION 
ON  THE 
TANK  WALL 

TRIAL 

NUMBER 

FOR  GAGE 

HEAT  LOSS 
PER  AREA 
BTU/FT  A  2 

WALL  AREA 
REPRESENTED 
IN*2 

TOTAL 
HEAT  LOST 
BTU 

1 

ft* 

0.2655275 

28.274 

0.0521356 

2 

E2 

0 . 4253386 

47.124 

0.1391921 

3 

E3 

0.4508400 

56.941 

0.1782728 

4 

D3 

0.3764253 

78.540 

0.2053086 

5 

D4 

0.2168532 

94.248 

0.1419304 

6 

D5 

0 . 3276899 

94 . 248 

0.2144730 

7 

C3 

0 . 4223368 

94.248 

0.2764194 

8 

C4 

0.2916098 

94.248 

0.1908586 

9 

E4 

0.3129244 

94.248 

0.2048090 

10 

B3 

0.4977912 

94.248 

0.3258043 

11 

B4 

0.3131344 

94.248 

0.2049465 

12 

B5 

0.4033936 

94.248 

0.2640211 

13 

F3 

0 . 4254639 

94.248 

0.2784661 

14 

F4 

0 . 2767622 

78.540 

0.1509507 

15 

F5 

0.3774505 

62.832 

0.1646942 

16 

F2 

0.3742415 

34.361 

0.0893008 

17 

B2 

0.3950619 

28.274 

0.0775693 

18 

C2 

0.3618079 

14.137 

0.0355200 

19 

D2 

0.3845312 

1.767 

0.0047185 

TOTALS  . . . 

.  1279.022 

3. 1993912 

ft 

For 

locations,  see  : 

ft* 

B1 

0.2407778 

Cl 

0 . 2777149 

D1 

0.2814635 

El 

0.2444963 

FI 

0.2831850 

page  50 . 


AVERAGE 


0.2655275 


Gage  left  in  same  spot  on 
tank  for  5  runs  to  verify 
experiment  repeatability. 
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